TRIPLE CORRELATION OF THE RIEMANN ZEROS 
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Abstract. We use the conjecture of Conrey, Farmer and Zirnbauer for averages of ratios of the 
Riemann zeta function [11] to calculate all the lower order terms of the triple correlation function 
of the Riemann zeros. A previous approach was suggested by Bogomolny and Keating [6] taking 
inspiration from semi-classical methods. At that point they did not write out the answer explicitly, 
so we do that here, illustrating that by our method all the lower order terms down to the constant 
can be calculated rigourously if one assumes the ratios conjecture of Conrey, Farmer and Zirnbauer. 
Bogomolny and Keating [3] returned to their previous results simultaneously with this current work, 
and have written out the full expression. The result presented in this paper agrees precisely with 
their formula, as well as with our numerical computations, which we include here. 

We also include an alternate proof of the triple correlation of eigenvalues from random U(N) 
matrices which follows a nearly identical method to that for the Riemann zeros, but is based on 
the theorem for averages of ratios of characteristic polynomials |12l I13| . 
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1. Introduction 



In 1973 Montgomery [22J proved the following result, assuming the Riemann Hypothesis, on the 
two-point correlation of the zeros of the Riemann zeta function: 
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for suitably decaying functions / with Fourier transform supported in [—1,1] and weight w{x) 
4 _ ) 4 m 2 ■ He conjectured that (|l.lj) would in fact hold for any test function /. 

In 1994 Hejal [16] proved a similar result for the triple correlation of Riemann zeros: 
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and the Fourier transform of the continuous, 



with weight w(x 1 ,x 2 ,x 3 ) = Y\ j<k exp [ - \{xj - x k ) 
suitably decaying test function / is supported on the hexagon with successive vertices (1,0), (0, 1), 
(-1, 1), (-1, 0), (0, -1) and (1, -1). Here S{x) = 

This was extended to the n-point correlation function and to more general L-functions by Rudnick 
and Sarnak [23] in 1996. 



These results encompass the rigorous work on the subject, but are limited by two things. Firstly, 
the support of the Fourier transform of the test function is always confined to a restricted range. 
Secondly, only the asymptotic for large T is found. This second point is understandable, as the goal 
was to show that this limiting form was the same as that for the n-point correlation function of 
eigenvalues from large-dimensional matrices from the GUE ensemble of random matrix theory (see 
[S] or |20] for review articles on the connection between random matrix theory and number theory) . 
This aim was duly achieved, but there is clearly interest in the lower-order terms, as Bogomolny and 
Keating's early results [B] showed that in the two-point correlation function of the Riemann zeros 
one sees sensitivity to the positions of the low Riemann zeros themselves - something that clearly 
does not happen in random matrix theory at any order. After it was predicted by Bogomolny 
and Keating, a numerical illustration of this for the two-point correlation function was first shown 
in [2], where Berry and Keating also fully explain a similar phenomenon in the number variance 
statistic first observed by Berry in 1988 [lj. A numerical plot of the two-point correlation function 
calculated using the first 100 000 zeros of the Riemann zeta function is shown in Figure [T] The 
x-axis is the correlation distance between pairs of zeros in unsealed units, showing the distinctive 
dip at each Riemann zero. The plot is a histogram of the number of pairs of zeros with a given 
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Figure 1. The two-point correlation function of the Riemann zeta function cal- 
culated from the first 100 000 zeros. The correlation distance is plotted along the 
x-axis. Note the significant dip at the first few Riemann zeros: 14.13, 21.02, 25.01, 



separation distance, and the y-axis is divided by T f J > where T ~ 75000 is the height of the 
100 000th zero. 

Meanwhile, alongside this rigorous work in the number theory community, physicists using semi- 
classical techniques applied in the field of quantum chaos, treated the Riemann zeta function as a 
model system (the prime numbers playing the role of periodic orbits) and so studied the correlation 
of Riemann zeros in analogous ways to those in which they would study correlation of energy levels 
in their more standard physical systems. The first step in this direction was the derivation by Keat- 
ing [T7] of the limiting form of the two-point correlation function of the Riemann zeros. This is a 
heuristic calculation and it relies on a conjecture by Hardy and Littlewood [15] for the behaviour of 
correlations between prime numbers, but it has the advantage that there are no restrictions on the 
support of a test function. Using analogous methods this result was then extended in two papers 
by Bogomolny and Keating [7] to obtain the limiting form of the n-point correlation function. 

The first result on lower-order terms of the correlations of the Riemann zeros was also by Bo- 
gomolny and Keating [6], but using a different heuristic inspired more directly from semiclassical 
methods. Here they truncate the Euler product for the Riemann zeta function at primes less than 
log (in semiclassical language this means taking periodic orbits up to the Heisenberg time) and 
define a new set of zeros from the resulting approximation to the staircase function of the Riemann 
zeros (the function that increases by one at the position of each zero). It is the two-point corre- 
lation of this new set of zeros that miraculously gives all the significant lower-order terms of the 
analogous statistic for the Riemann zeros. Using the same method, an expression is also given [6] 
in semiclassical language for all lower-order terms of the three-point correlation function, and it 
is this which could be turned into a detailed formula including all terms calculated in the present 
paper, but the authors did not publish it explicitly at that time. 
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Bogomolny and Keating also did further work, largely unpublished, obtaining all the lower-order 
terms for the 2-, 3- and 4-point correlations both by extending their Hardy-Littlewood method (see 
[18] for some details of the two-point correlation function) and by the method mentioned in the 
previous paragraph, as well as two other heuristic methods (see [3]). 

Recently [3] they have written out the lower-order terms for the three point correlation function 
of the Riemann zeros in full detail and these agree with the results presented in this paper. 



2. The Riemann zeros 

2.1. Results. Using the ratios conjecture of Conrey, Farmer and Zirnbauer \12\ 111] we obtain the 
following 

Theorem 2.1. Assuming the ratios conjecture and summing over distinct zeros of the Riemann 
zeta function: 
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The expression (|2.ip is plotted in Figure [21 with f(v±,V2) = 8{v\ — x)5(v2 — y)- The x- and y- 

/ log \ 3 

axes are unsealed, but (|2.ip is divided by T( 2 ^ J . The density plot is light for large values and 
dark where (|2.ip is small. Note the horizontal, vertical and diagonal lines occurring at the Riemann 
zeros, caused by terms like ^-(1 + ix), ^-(1 + iy — ix), etc. Numerical computation of (|2.ip breaks 
down near the x- and y- axes and on the diagonal because of the principal value integration, so 
the plot has been set to zero in these regions. The plot could be completed with a more careful 
expansion of the formula around x = 0, y = and x = y, but this would not be particularly 
edifying. The result would be extremely similar to the random matrix limit shown in Figure 
The maximum height on the contour plot in Figure [2] is about 0.799. An idea of the height of 
the plot can be seen in Figure [31 which is a horizontal cross-section of Figure [2] at height 5 on the 
y-axis. 

In Figure [4] the numerical triple correlation, using the first 100 000 zeros, is plotted, again scaled 
as above. The difference between this and Figure [2] is shown in Figure 02 where the maximum 
height of the plot is about 0.164. Compare this with the maximum height of Figure HJ which is 
about 0.923. The maximum value of the difference plot may seem rather large, as we expect an 
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Figure 2. The triple correlation of the Riemann zeros: we plot the quantity in 
square brackets from (|2.ip . The x and y axes are v% and V2 and the density plot is 
lighter were the function has a higher value and is darker for smaller values. Note 
the horizontal, vertical and diagonal lines occurring at the Riemann zeros: 14.13, 
21.02, 25.01, ... 



error of T 1 / 2+€ , but this is probably due to the relatively small value of T used for these plots. 
For T ~ 75000, T~ 1 / 2 is around 0.003, but for T values of this size powers of logT can make a 
big difference. The mean value of the points on the difference plot (Figure [5]) is -0.00127, and the 
standard deviation is about 0.03, which gives a better idea of the spread of the points. We also 
note that the mean of the absolute value of Figure [5] is about 0.0257. 

For ease of comparison with Theorem 12. 11 we state here the similar random matrix result which 
is derived in detail in Section 13.21 The identical structure of Theorem 12.11 and Theorem 12.21 is 
apparent if we recall the equivalence N = log (see for example [19] ) and reduce (|2.8p from three 
to two variables by performing a simple translation such as 02 — > 02 + 0\ and #3 — > #3 + 0\ and 
noting from fl2J)D and (pTEO]) that J(a,f3;j) = J(a + t, /3 + t;j -t) and J{a;f3) = J(a + t; (3 -t). In 
comparing J[a.,f3;^) with I(a,fi; 7) and J{a;(3) with Ii(a;/3), z(x) = (1 — e~ x )~ l plays the part 
of C(l + x), as is always the case in moment and ratios conjectures. Using the ratios theorem of 
Conrey, Farmer and Zirnbauer \\.2\ [TTj we obtain the following 
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Figure 3. A horizontal cross-section of Figure [2] at height 5 on the y-axis to illus- 
trate the height of that plot. 



Theorem 2.2. With the star indicating a sum over distinct eigenvalues we have: 

W)--= [ Y, f(o r .o n .o !: ,uix 



where 



(2.9) 



1 



(27T) 




J(i9±,i9 2 ;-i9 3 ) 



+J{i9 1 ,i9 3 ; -%9 2 ) + J{i9 2 , i9 3 ; -i9 x ) + J(-i9 h -i9 2 ; i9 3 ) 
+J(-i9 1 , -i9 3 ; i9 2 ) + J(-i9 2 , -i9 3 ;i9{) 

+N(J(-i9r,i9 3 ) + J(-i9 2 ; i9 3 ) + J(-i9 1 ;i9 2 ) 
+J(-i9 3 ;i9 2 ) + J(-i9 2 ; i9 x ) + J(-i9 3 ;i9 1 )) 

(2.8) +N'Af(9 1 ,9 2 ,9 3 ) d9 x d9 2 d9 3 



J(ai,a 2 ;f3) 



e 



Ju(N) Ax Ax Ax* 

-N(<*i+P) z ( ai + f3)z(-ai - (3) (^(a 2 - ax) -^(a 2 + 0)\ 



8 



J.B. CONREY AND N.C. SNAITH 




FIGURE 4. The numerically calculated triple correlation of the Riemann zeros. The 
distribution of triplets of zeros is depicted by plotting the distance between the first 
and second zero in the triplet on the horizontal axis against the distance between 
the first and third zero on the vertical. Higher occurrences appear lighter grey and 
where no triplets fall the plot is black. 



and 

J(a;/J): = e"W ^'^{e^dX 
JU(N) A X Ax* 

(2.10) = (^\ (a + P) + e- N ^z(a + P)z(-a-P). 

2.2. Moments of the logarithmic derivative of the Riemann zeta function. Calculating 
correlation functions of the Riemann zeros using the conjectural formulae for averages of ratios of 
zeta functions proceeds via moments of the logarithmic derivative of the Riemann zeta function: 
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Figure 5 . The difference between Figures [4] and [2j 
A property of these moments that will be of use to us later is that 

(2.12) M a i> • • • > a k;/3i, ■ ■ ■ = iriPx, ■ ■ ■ ,0e;®l, ■ ■ ■ ,®k)- 

For now the arguments of I r will always have positive real parts. In this case we have an 
approximate translation invariance: if ^(aj), ^R{(3j), 9?(aj + A), — A) > then 

(2.13) I r (ai + A, . . . , a h + A; ft - A, . . . , ft - A) 

= / r (ai,...,a fc ;ft,...,ft) + 0(|A|T e ). 

This can be seen by a change of variables t —> t — iX in (|2.1ip and using RH to bound C'/C( c7 + 
by t e . 

In particular, to calculate the three-point correlation function, we will need formula (12. 2ft for 
/(ai, a2; /?) (here we introduce the convention that / = Iq). 

We will now proceed to derive that formula using the form of the ratios conjecture |llj for three 
zeta functions over three zeta functions (with the conditions — \ < 3lai, 3iai2, < \, ^ 
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3f?7i, 318 < j and So^, 3/3, Qtyi, 372, 9tf <^ e T 1_e ): 
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The main term of (|2,14p is analytic in the specified range of parameters; the apparent poles from 
Y cancel, as can be checked directly, or by writing the three terms as a contour integral as in [11] 
or Section 2.5 of [TO] . 

Here fJ,(n) is the Mobius function and the final expression above reflects the fact that fi(p m ) is 
1 for m = 0, it is —1 for m = 1 and zero for any power of p higher than the first. Also note that 
yl c (ai,a 2 ,/5;7i,72,/3) = 1- 

We will not go through the reasoning behind this conjecture in full, but a simpler example (one 
zeta function over one zeta function) can be seen in full detail in [14J and the original reference for 
the general case is [11]. We will say only that the recipe for arriving at a ratios conjecture involves 
replacing each zeta function in the numerator with an approximate functional equation, those in 
the denominator with the Dirichlet series for and then applying a series of rules to discard all 
the terms in the resulting multiple sums except for those seen above in At. The purpose of Y is to 
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factor out the divergent terms in the these sums, leaving Aq convergent for small values of a%, a 2 , 
/3,7i,72 and 5. 

Armed with this ratios conjecture, we want to evaluate 
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Examining the derivative of the first term in (|2.14j) . we find a great deal of cancellation upon 
substituting 71 = «i, 72 = a 2 and 6 = (3, and the only surviving terms are 
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performed on these factors. The final differentiation with respect to a 2 shows us that (j2.20j) equals 
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Substituting these expressions into (12.190 . (I2.20p and (I2.23D . we arrive at (|2.2p as expected. 

The other version of (|2,lip that we need is I\{a;(5), as given in (|2.3p . This calculation has in 
fact already been carried out in |14| using the two zeta functions over two zeta functions ratios 
conjecture in a manner very similar to the three-over-three calculation above. That is |X2|, 111], with 

-\ < 3fta,3fy3 < |, < 3^7,^(5 < \ and 9a, 9/3, 97, 95 <C e T 1-6 (for every e > 0), the ratios 
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provided that <C 9?a,9?/3 < |. The ratios conjecture recipe can incorporate the log ^ factor 
in Ii(a; (3) without any alteration, giving (|2.3p . 



2.3. Triple correlation as a contour integral. We start with the a triple sum over zeros of 
the Riemann zeta function, where the zeros do not have to be distinct. Using Cauchy's residue 
theorem, the triple sum can be written as a triple contour integral where each contour is a rectangle 
enclosing the zeros 1/2 + 17 (assuming the Riemann Hypothesis) with heights < 7 < T: 

(2.31) 5(71,72,73) 

0<7i.72,73<T' 

jg(-i{x - 1/2), -i(y - 1/2), -i(z - l/2))^(x)^(y)^(z)dxdydz. 




In investigating the 3-point correlation of the Riemann zeros, we are interested in the relative 
spacing between triples of zeros, so we can assume that the test function is translation invariant. 
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Thus we define a function / satisfying the conditions 

(2.32) f(x,y) is holomorphic for < 2, \Qy\ < 2 

and that f(x,y) < 1/(1 + \x\ 2 + \y\ 2 ) as \x\ or \y\ — > oo. 

Thus we have 

/(7i "72,71 -73) = I I I '^(zi)^-(z 2 )^(z 3 ) 

0<7i, 72,73 <T 

(2.33) x/( - i(zi - z 2 ), -i^i ~ Z3))dzidz 2 dz 3 , 

where the contours are rectangles with corners at the points (a, 0), (a, iT), (1 — b, iT) and (1 — b, 0) 
with 1/2 < a, b < 1. We distinguish a and 6 only for ease of following the manipulations in the 
following calculations. 

The horizontal portions of the contour of integration can be chosen so that the integral along 
them is negligible, so we concentrate on the vertical sides of the contours. This makes (|2,ip the 
sum of eight integrals J\, . . . ,Jg which will be met one by one below (a subscript a on an integral 
indicates integration from (a, 0) to (a,iT)). 



■h -r-T / / / T(*l)Tte)T(«0 



First we have 

(27H) 3 JaJaJaC^ C 

x/( - i{zi - z 2 ), -%{z\ - z 3 ))dz 1 dz 2 dz 3 
(2.34) = 0(T € ). 

The final line is true because all three contours can be moved to the right (assuming the Riemann 
Hypothesis) where ^- converges and can be integrated term by term (the pole at 1 doesn't contribute 
more than a constant). 



h ■ -T— T7 / / / ^(Zl)^(z 2 )^(z 3 ) 



Next we examine 



(2vri)3 JaJJ^C LJ C" J C 

(2.35) x/( - i(zi - z 2 ), -i(zi - z 3 ))dz x dz 2 dz 3 . 

We use the functional equation 

(' x' (' 

(2.36) L (s) = K (s) _L(i_ s ) 



and obtain 



C w x w C 
1 f f f fx' , ,W, X. 



j2 = -^yjj M (|w-^-^)W w 

(2.37) x/( - i(zi - z 2 ), - z 3 ))dzxdz 2 dz 3 . 

For the term with ^r(zi) we can shift the integrals to the right, as we did with Jj, and so see that 
the contribution is only 0(T £ ). A similar manoeuver cannot be done with the integral containing 
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the three however, because the z\ integral is to the left of the critical line. So, 
1 f T f T f T (' (' (' 

j2 = wrL I I c (6 -^c (a+rf2) c (a+ ^ 3) 

(2.38) x/(ti - t 2 - i(l - b - a),ti - t 3 -i(l-b- a)) dhdt 2 dt 3 + 0{T € ). 
Now, letting v\ = t 2 — h and v 2 = t% — 1\, we have 

xf{-v 1 -i(l-b-a), -v 2 — i(l — 6 — a^dv^dh + 0{T e ) 
= (27F J_ T J_ T f(-v 1 -i(l-b-a),-v 2 -i(l-b-a)) 

(2.39) x / y(6-i*i)y(a + i(ui + ti))~(a + ifa + ti))dtid«idv 2 + 0(T e ). 

In this last line we have switched the order of integration of ti with v\ and t>2- The range of the 
innermost integral should really be from max{0, — v i, — v 2 } to min{T, T — Vi 3 T — v 2 }. However, 
since we are assuming that f(x,y) decays fast, see (|2.32[) . it is not hard to show that extending 
the interval to (0, T) as we have done above, and will do again in the following integrals J3 to J$, 
incurs only an error that is a power of logT and so can be incorporated into the error term 0(T € ). 

Now all that is left is to tidy up the expression for J 2 , so we shift the contours of integration off 
the real axis and apply the definition of I(a±, a 2 ; (3) from (|2.11j) with r = 0: 

j pT-i(l-b-a) pT-i(l-b-a) 

J 2 = 7^3 / / f(.Vl,V 2 ) 

{Zir)° J- T -i(l-b-a) J-T~i(l-b-a) 



X 



r T c c c 

/ -y(b- ih)jr{l + ih-b- +ih-b- iv 2 )dt 1 dv 1 dv 2 + 0{T e ) 

Jo C C C 



l 



T-t(l-6-o) r T-i(l-b-a) 



(2.40) = / / f(v 1 ,v 2 )I(-iv 1 ,-iv 2 ;0)dv 1 dv 2 + O(T e ). 

(27Tj° J~T-i(l-b-a) J-T-i(l-b-a) 

To simplify the last line we have used (|2.13p . 

We use exactly the same sequence of manipulations to obtain similar expressions for J3 and J4: 



(2my J a J a c C C 



x/( - ~ z 2 ), -i(zi - z-i))dz x dz 2 dzi 

pT pT+i{l-b-a) 

(2.41) = / / /(v 1 ,« 2 )J(0,-i« 2 ;iT; 1 )dt;idt;2 + O(r £ 

-T J -T+i(l-b-a) 



(2tt) 2 
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and 



(2mf A_ 6 JaJaC T '( 

x/( - i(zi - z 2 ), -i{z\ - z 3 ))dz 1 dz 2 dz 3 

i rT+iil-b-a) r-T 

(2.42) = / / f(v 1 ,v 2 )I(0,-iv 1 ;iv 2 )dv 1 dv 2 + O(T e ). 



The integral J5 throws up something slightly different. We start off in an identical manner, 
replacing ^-(^2) and ^-(z 3 ) by their functional equation (|2.36j) . This results in four terms, one of 
which contributes no more than 0(T € ). We are left with three terms: 

^5 := tAt III 



(2^) 3 Jl-bJl-bJa C C C 

x/( - i{z\ ~ z 2 ), -i(zi - z 3 ))dz 1 dz 2 dz 3 

1 r T r T r T ( c c c 

-— (1 - b + 1*2)77(0 + ihA(b - it 3 ) 

x C C 

- — (1 " b + tt 3 )£(a + iti)£(b - it 2 )) 

x C C J 

x/(*l - t 2 + i(l-a- b),ti -* 3 + i(l - a - b))dtidt 2 dt 3 + 0(T e 

1 /" T f T ~ tl f T ~ h ( (' c c 



(2tt) 



y' C C 

- — (1 - 6 + i( Vl + *i))-7(a + i*i)-7(6 - i(v 2 + *i)) 

X C C 

--(1 - b + i(v 2 + h))^-(a + iti)-V(6 - + *i)) J 

x c c y 

(2.43) x/( - vi + i(l - a - 6), -u 2 + i(l - a - b))dv!dv 2 dti + 0{T e ), 

where in the last line we have made our usual change of variables v\ = t 2 — *i and v 2 = t 3 — t\ 
We now note the asymptotic for — : 

(2.44) ^( 1/2 + it) = _l„ g M( 1 + o(i [ 
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Since / is very small when v\ or v 2 are large, we replace — (1 — b + i(vj + z{)) with — log and 
obtain 



/•T rT 

(2vr)=" ' 



J 5 = To~^3 I t J f[-vi + i(l-a-b),-v 2 +i(l-a-b)) 

C C 
M6-»(vi+*i))^r(6-i(«2 + ti)) 



+ log |^(a + iti)^(6 - i(«2 + *i)) 



+ log |^(a + - *(«i + ti)) ] dhdvxdvi + 0(T £ ), 



where we have extended the range of ii, after exchanging the order of integration, by the same 
argument as for J 2 . 

Now we can write J5 as a double integral along contours running just below the real axis, in a 
similar form to J 2 , J3 and J4 above, but we introduce Ii(a;/3) (see (|2.1ip and (|2.3p ) to arrive at 

^ rT+i(l-a-b) i-T+i(l-a-b) 

J 5 = TTTTT / / /(«1,U 2 )( J(0;Wi,TO 2 ) 

(^7TJ J 1 /_ T+i ( 1 _ a _ 6 ) J-T+i(l-a-fe) V 

(2.45) +Ji(0; tua) + li(0; wi)Wdu 2 + 0(T £ ). 



Proceeding exactly as for J5, we obtain similar expressions for J§ and J7: 



x/( - i(zi - z 2 ), -i{z\ - z 3 ))dzidz 2 dz 3 



1 



T r T-i(l-a-b) 

f(vi,v 2 )(l(-ivi;0,iv 2 ) 



(2-7T) 3 J_ t J-T-i(l~a-b) 

(2.46) +/i (-ivi ; iu 2 ) + 7i(-iui;0) )duidu 2 + 0(T e ) 



and 



(2vri)3 JJ 1 _J 1 _ b ( K J C 

x/( - «(^l - *2), -iC^i ~ Z3))dzidz 2 dz 3 

I rT-i(l-a~b) rT 

= 7o~TT / / /(^^(/(-w^O,^!) 

(2vr)3 y_ T _ i(1 _ _ 6) y_ T V 

(2.47) +/i(-*V2;»«i) +/i(-iu 2 ;0))di;idu2 + 0(T e ). 
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We are just left with the integral 



^8 := -tA^ / / / ^1)^2)^3) 



(2ni) 3 A-& A-& A-6 C C C 

(2.48) x/( - i(zi - z 2 ), -i(zi - z 3 )) dz 1 dz 2 dz 3 

to evaluate. Once each ^(z) is replaced by its functional equation (|2.36p . any term with at least 

one ^-(1 — z) in it can be shown to be size just 0(T e ) by shifting the contour far to the left. This 
leaves us with 

(2.49) x/(ti - t 2 , h - h)dhdt 2 dh + 0{T e ). 

Here the contours have all been moved onto the half-line as there are no longer zeros of zeta 
functions to avoid. We replace t 2 — t± with v\ and replace t 3 — t\ with v 2 and substitute — log |^ 

for each 2L factor with the help of (|2.44|) . as we did in the discussion of J5, and so obtain 

A 

(2.50) J 8 = J T ^J T J ( Vl ,v 2 ) J ^ log* ±dtd Vl dv 2 + 0(T% 



The sum of integrals J\, . . . , Jg gives us (see ()2.53[) ) the expression for 
(2-51) ^2 /(7i -72,71 -73) ■ 

0<7i.72,73<T 

Such a triple sum over zeros necessarily contains terms where two zeros are identical. These are 
essentially two-point correlations, rather than three-point statistics. Similarly, terms where all three 
zeros are identical are just one-point statistics. To remove these lower-order correlations, one looks 
at the sum 

(2.52) Y, /(7i "72,71 -73) 

0<7i^72^73<T 

instead. As will be shown in the next section, simply rewriting as an integral on the real line 
each of the integrals above that occurs with a contour on (— T — i(l — a — b),T — i(l — a — b)) or 
(— T + i(l — a — b),T + i(l — a — b)) gives us the purely triple correlation, and so we obtain the final 
result (12. ID . The contributions from the poles that we meet as the contours are shifted to the real 
axes are the source of the two- and one-point correlation terms. 
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2.4. Contributions from lower-order correlations. Collecting together the results from the 
previous section we have 

/(7i -72,7i -73) = 7^73 f J J f(vi,v 2 ) J log 3 ±dt dv x dv 2 

/T+i(l-o-6) /•T+i(l-a-6) 
/ /(wi,u 2 )( J(0;iUi,iU2) 



T-H(l-o-6) J-T+i(l-o-6) 

+Ji(0;ro 2 ) + -fi(0;zui) )eh;ich; 2 



/ /(«i,w 2 )( J(-wi;0,w 2 ) 

-T J —T—i(l—a—b) V 

+Ji(— iui;iu 2 ) + ii(— ivx;0)jdvidv 2 

/T-i(l-a-fe) /-T 
/ /(«i,w 2 )( J(-w 2 ;0,wi) 
-T-i(l-o-b) J-T v 

+Ji(— m; 2 ;u>i) + Ji(— «u 2 ;0)Jdt»idu 2 

«T— i(l-6-a) pT-i(l-b-a) 

I f(vi,v 2 )I(-ivi,-iv 2 ;0)dvidv2 
' -T-i{l-b-a) J-T-i{X-b-a) 

/T rT+i(l-b-a) 
/ f(vi,v 2 )I(0,-iv 2 ;ivi)dvidv2 
-T J -T+i(l-b-a) 

fT+i(X-b-a) f-T \ 

(2.53) +/ / f(vi,v 2 )I(0,-iv 1 ;iv 2 )dv 1 dv 2 ) +0(T € ). 

J -T+i(l-b-a) J-T J 

The goal of this section is to move all the contours of integration above onto the real axis, 
evaluating the contributions from the poles encountered during this process. The resulting integrals 
along the real axis will then be computable, for the purposes of Figure ???? for example, as principal 
value integrals. Elegantly, the terms resulting from the residue at the various poles will yield the 
contribution to (|2.53p from lower-order correlations between Riemann zeros. 

Consider, as an example of the method, the term (coming from J5) 

/T+i(l-a-b) rT+i(l-a-b) 
/ f(vi,v 2 )l(ivi,iv 2 ;0)dvidv 2 
-T+i(l-a-b) J -T+i(l-a-b) 

from (|2.53j) . (Here we have switched the order of the arguments of / using (|2.12j) .) Treating the 
integrand first as a function of v 2 alone (leaving v\ fixed and non-zero), we see that the v 2 contour 
lies below the real axis (since a > 1/2 and b > 1/2) and there is a pole at v 2 = 0, so shifting the 
contour onto the axis results in a principal value integral in v 2 plus itr times the residue of that 
pole. Using the definition of the sums and products over primes found in (12. 4ft to (12. 7j) . expanding 
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around v 2 = gives us 

(2.55) Q(ivi,iv 2 ) = O(l) 

(2.56) A(iv 2 ) = l + 0(t|) 

( 2 -57) P(iv2,m) = iv 2 ^ 1)2 + 0(^ 2 2 ) 

(2.58) P(iui,w 2 ) = O(l). 

These, in conjunction with the poles of l(iv\,iv 2 ;0) at v 2 = resulting from ("(1 — iu2)C(l + iv 2 ) 
and ^-(1 + iv 2 ), allows us to express (12.541) as 

fT rT+i(l-a-b) 



f(vi,v 2 )I(ivi,iv 2 ; 0)dv\dv 2 

-T J -T+i(l-a-b) 



pT+i(l-a-b) pT / a/. , + 

+7T / ( (y\ (1 + ivi) + (^-) C(l + wi)C(l - ivx)A{w x ) 

J-T+i(l-a-b)J0 VH 7 27T 

( 2 - 59 ) " E (pi-Sf-iy. ) 0)* **■ 

where the t>2 integral is understood as a principal value integral. 

Now we move the v\ integral onto the real axis. In this case also we encounter a pole at V\ = 0. 
The first integral in (|2.59j) yields a residue contribution essentially identical to that from the v 2 
pole treated above, and the second integral has a pole at v\ = with residue — ilog We take 
in times the contributions from these poles and the end result is 

/ / f(vi,v 2 )l(ivi,iv 2 ;0)dvidv 2 
J-T J -T 

+7T J T f( v i,0)f ((^)\l + ivi) + ^y W \(l + ivi)£(l-iv 1 )A(iv 1 ) 

-E (p i+^!i )2 )^^ 

+^ y ' /(0,«2) ^ (Y^)'(l + w 2 ) + (^)^ 2 C(l + ^ 2 )C(l- i«2)A(i«2) 

-E<^£! 1)a )*** 

(2.60) +ir 2 [ T f (0,0) log ±dt. 

Jo 

We note that the final term comprises information on the one-point correlation function, whereas the 
second and third integrals in (|2.60p will form a contribution to the two-point correlation function, 
as will be discussed at the end of this section. 



TRIPLE CORRELATION OF THE RIEMANN ZEROS 



21 



In the meantime we will perform sample calculations on three other terms from (12.530 . leaving 
the remainder of the terms to the reader, since they all follow a similar pattern. 

We consider now 

/T+i(l-a-b) i-T+i(l-a-b) 
/ f(vi,v 2 )h(0;iv 2 )dv 1 dv 2 . 
-T+t(l-o-6) J-T+i(l-a-b) 

Examining expression (|2.3H for li(0;iv 2 ), we see that the (^)' term has a second order pole at 
V2 = with residue zero, the term containing £(1 + zu2)£(l — iv 2 ) has residue 

T 



(2.62) / -Hog 2 ^dt, 

Jo 



and B(0,v 2 ) is analytic near = 0. The integrand is analytic in v±, so moving that contour onto 
the real axis does not incur any polar contribution. Thus (|2.6ip equals 

(2.63) f [ f(v 1 ,v 2 )I 1 (0;v 2 )dv 1 dv 2 + TT [ /(«i,0) / log 2 ±dt dv x . 

J-TJ-T J -T Jo 

The next term that deserves consideration is 

/T-i(l-a-b) pT 
\ f(vi,v 2 )l(ivi,0;-iv 2 )dv 1 dv 2 , 
-T-i(l-a-b) J-T 

where again we have used (|2.12p to exchange the arguments of I preceding the semicolon with those 
following it. We note that in these integrals where two consecutive principal value calculations 
are made the order of the integrals cannot be exchanged after the first principal value integral 
is obtained, so we always address the principal values starting with the outermost integral and 
working inwards. 

An inspection of (|2.2p reveals that I(ivi,Q; —iv 2 ) is not singular at v\ = (the poles of two ^ 
terms cancel), but as we move the v 2 contour onto the real axis we encounter extra difficulties when 
V\ = 0, so we will start by temporarily shifting the v\ contour so that it runs just below the real 
axis; we choose below rather than above the axis so as to avoid the pole at v\ = v 2 . 

Now that the v\ contour does not pass through zero, we can move the v 2 contour onto the real 
axis and pick up exactly the same polar contribution as when evaluating (|2,54p , (An extra minus 
sign in the residue compensates for the fact that this time we need it multiplied by —in since we 
are half-circling the pole in the clockwise direction due to the original contour passing above the 
real axis.) Thus, for some e > 0, (|2.64p is 

(■T f-T-ie 

f(vi,v 2 )l(ivi,0; -iv 2 )dv\dv 2 

T / '('\' ( t \ 



-T J -T-ie 



+7T / / ^ )(l+i Vl )+ — ({1 + iVl)((l - i Vl )A(iVl) 



T-ie Jo \^C J V2 7T 

( 2 -65) - V j-M^ ) /(«!, Q)dt dm, 



where the v 2 integral is understood as a principal value integral. 
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As we move the v\ contour back to the real axis, we encounter a pole of the first integral above 
at v\ = v 2 (the apparent pole at v\ = in fact cancels), and a pole of the second integral at v\ = 
with residue — ilog-^, as before. The final result is that (|2.64|) equals 



T r-T 

\ f(vi,v 2 )I{iv 1 ,0;-iv 2 )dvidv2 

-T J-T 

+vr f( V1 , o) £ ((^)'c 1 + ^i) + (^r™ 1 ^ 1 + ivi ^ - ^imk) 



-E (p i-S! 1)2 )^ 

(2.66) +7T 2 F f (0,0) log ±dt. 

Jo 



The last term we will consider (all the others follow in an identical manner to one of those 
discussed) is 



/T i-T-i{l-a-b) 
\ f(vi,v 2 )l(iv 2 ,0;-ivi)dvidv 2 , 
-T J -T-i(X-a-b) 



where as usual we have used fj2. 12f) to exchange the arguments of / preceding the semicolon with 
those following it. We write this as 



/[■T-i(l-a-b) 
\ f(v±,v 2 )l(iv 2 ,0;-ivi)dvidv 2 . 
-SVT] J _ T _ i(1 _ a M 
\v 2 \>e V I 



e- 



This has not changed a thing because the outer integral is perfectly well-behaved at v 2 = 0, but it 
means that as we move the inner integral onto the real axis we avoid the complications that arise 
if v 2 = 0. 
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The inner integral has poles at v\ = and v\ = i>2 with the usual and now-familiar residues, so 
(gHH) equals 

rT r T 

f(vi,V2)I(iv2,0; —iv\)dvidv2 



T J — T 

T 



+7T J J(0,V 2 ) J ((^)'( 1 +^2)+ (^) m \{l + iv 2 )C{l-iV2)A{iv 2 ) 

-E ( pi+S! 1)2 )^^ 

+7T j^f(v2,V 2 )j^ (^(^y(l-iV2)+^y V2 ((l+iV2)((l-iV2)A(-iV2) 

(2.69) -V ; : lQg2p , )dtdv 2 , 

where the integrals in v% and t>2 are to be interpreted as principal value integrals. 

All of the other terms in (|2.53p can be handled in exactly the same way as the one of the four 
treated here. The complete result is 

rT r T 



0<7ii72,73<T 



E / (Ti - 72, 71 - 73) = ^3 J t f t f(vi,v 2 ) 

I / log 3 —du + I(ivi,iv2; 0) + 1(0, iv\; -iv 2 ) + 1(0, iv 2 ; —iv\) 

-if 2 ; 0) + /(0, — iv2] ivi) + 1(0, —ivy, iv 2 ) 
+li(0;iv 2 ) + h(0;ivi) +/i(-w 2 ;wi) + h(-iv 2 ;0) + I\(—ivi; iv 2 ) 



+li(—ivi;0) \ dv\dv 



(2^y 



2 



T rT rl 

f(0,v 2 )s(v 2 )dv2 + / f(v 1 ,0)s(vi)dv 1 + f(vi,vi)s(-vi)dvi 

J-T J-T 



T 



rT f-T 

+ I f(0,v 2 )s(-v 2 )dv2 + / f(vi,0)s(-vi)dv 1 + I f(vi,vi)s(vi)dvi 



(2-70) /(0,0) log £di, 



where 



f T 1 n /(T'x ' / ( \ -ix 

s(x) = / - log 2 £ + (i-J (1 + ix) + (— ) C(l + «)((! - 
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Note that some terms above can be combined if we include the natural assumption that f(x,y) = 
f(—x, —y). All the integrals on the interval (— T, T) should be considered as principal value integrals 
near any poles at the origin or at V\ = V2- 

Now we want to identify the terms in the last three lines of (|2,70p as lower order correlations. 
To do this, note that we can rewrite the triple sum over zeta zeros as 

f^ 1 ~ 72 ' 71 ~ 73) = ^ fin - 72, 71 - 73) 

0<71,72,73<T 0<7i 7^72 7^73 <T 

+ ^2 fin - 72, 0) + J2 f(°> 71 - 73) + ~ 72, 71 - 72) 

0<71^72<T 0<7it^73<T 0<7i^7 2 <T 

(2.72) + /(°>°)- 

o<7i<r 

The standard result on the density of the zeros of the Riemann zeta function gives, 

< 2 - 73 > E 1 -s"*sf- h + °«' 

0<7i<T 

for large T, and an expression for the two-point correlation, derived from the ratios conjecture, is 
given in 



<*■«> E /h - V) = <^ £ f / w ( * s + (7) ' (1 + 



77^7' <T 



+ f ^ C(l - *r)C(l + *r)A(»r) - B(tr) + ^ (1 - ir) 
+ ^^-j C(l - »r)C(l + *r)A(-»r) - B(-ir)^ (it dr + 0(T 1/2+e ), 

where /(z) is holomorphic throughout the strip \Ssz\ < 2, is real on the real line and satisfies 
f(x) < 1/(1 + x 2 ) as x — > oo. In addition, the integral is to be regarded as a principal value near 
r = 0, 

(2.75) AW = n (1 " 5 ^ )(1 "' + 5 ^ ) 



(i-i) 



and 

<™> b "»-Eto 



Op 1 

p 



We see immediately from (|2.73p . (|2.74p and (|2.72p that y^ 0<71 ^ 72 ^ 73<T /(7i ~ 72, 7i ~ 73) is 
given by the first five lines of (I2.70P , and this is the result presented in Theorem 12.11 
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2.5. Retrieving the asymptotic result. We want to confirm that our formula for the triple 
correlation of the Riemann zeros (|2.ip tends to the limit (plotted in Figure [6]) 



(2.77) 



hm -=- 



E 



I ('^(71 -72),^ (71 -73) 



2tt 0<7i^72^73<T 



2vr 



2vr 



OO /"OO 



oo J — OO 



OO /"OO 



OO J — OO 



f(vi,v 2 

i 

f(vi,v 2 ) [ 1 
+2 



1 5(ui) S-(«2) 

S(-ui) 1 S^fi-i^) 

S(v 2 ) S(vi - v 2 ) 1 



dv-\ dvo 



sin 2 (7r(fi — v-ij) 



sin 2 (7rui) 

2^ 



7T 2 (fl — V 2 ) 2 TT*V 

sin (7rui) sin (irv^ sin (7r(t>i — v 2 )) 

TTVi TTV 2 7r(vi — V 2 ) 



sin [TTV 2 J 
7T 2 Un 



| dv± dv 2 . 




Z.S 5. T.S 10. 12.5 13. 17 . 5 20. 



Figure 6. The triple correlation of eigenvalues of random matrices from U(N) in 
the limit as N —* oo. That is, we have plotted the 3x3 determinant from the second 
line of (|2.77p . Note that to compare this picture to Figures [2] and H] the axes of those 
figures would need to be scaled by the mean density of zeros, (log £-)/{2ir). 
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Using (|2.ip , with L = log ^ , we scale the variables in the test function by ^ , and make a change 

2nv 



of variables v — > =r in the vi and v 2 integrals: 

£ /(^(7i- 72), ^(71 -73)) 



TL TL 
1 f 2?r /" 2tt 



0<7i 7^72 7^73 <T 

"_L TL 

L 2->- I TL / ^Z^) 
2vr 2vr 

(2.78) xf [ T log 3 ^-du + I (^,^;0) +/(0,^L;-^p) + /(o,^p;-^i) 

\j( 27rtt)i 27TTO2 . r>\ 1 r /pi 27rif2 . 2iriv\ \ j_ T ( n 27rii>i . 2_7ra>2 "\ 

~ i ~ J I L ' L » u y V, ! L ' L / v 5 L 1 L ) 

+/i(0;^) +/i(0;^) ^p;0) 

+Jl ( _ 2zrp ; 2^2) + 7l ( _ 2^1; )^J + (T £ ), 

Starting with the first term of (|2.78p . we see that asymptotically for large T, 

TL TL 

1 f 2n f 2?r f 1 , i 

J _TL. j„TL f{Vl ' V2) J 2^**** 
t~i m /"CO /»oo 

(2.79) r "27r l0S 2vr' ' Si^i^dvidv^. 



00 ^ —00 



A little more work is needed for the other terms of (|2.78|) . Next we consider 

TL TL 

(2-80) -4" f(n,v 2 )l(^,^-,0)dv ld v 2 . 

J ~ 2tt J ~ 2tt 

Using the definition of I (ax, a 2 ; 0) in (|2.2|) . we will be a little imprecise and discard any terms that 
will not ultimately contribute to the leading-order Tlog^ result. We thus keep only terms with 
third-order poles as T — > 00, and the leading-order contribution to (|2.8U|) is contained in 



TL TL 
2n f 2n 



L 2 2tt ]_TL J_IL 

2tt 2tt 



f(vi,v 2 ) 



X 



I _ \ L /■ f-\ 2-Kivi \ /■ (-1 1 2nivi \ a ( 2iriv\ \ *5 ( -i , 2iriv2 2mv\\ 

V2vr/ ' Z v ^ L / V L J 7"^- L ~) 







] f J L 27TW2 ^^ I 2?r ™ 2 ) ^ 2-7TW2 ^2 (^ I 27Tt-»l 27TTO 2 ^ 



_ 2^ 2 , 

I __ I L ^/l _ 2-Ktvi \/-(-\ 1 27m)2 \ 4 / 27r»i) 2 \ 2_ 1 



(2.81) -( — ) L C(l-^)((l+ 2 ^)A(^p)^(l + ^i) )dtdv 1 dv 2 . 
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Keeping only the polar terms and performing the integral over T of the form 



(2.82) J g [-) - .»-,', 

we find (I2.80P is asymptotic to 



T , J. 2TTiv 



-L / /fa^ 2 ) [- — g— — + 



2vr J^ooJ-oo ' V *(2tt) 3 v 2 (u 2 - vi) i(27r) 3 u^ 2 



g- 2 7TO 2 g — 2 7ri«2 



(2.83) ~ 77Z ,n 2 , 7 + - fn ^ 2 dv l dv 2- 

A similar calculation can be done for 



TL 


TL 


2tt 


rw 


^ J 


1 TL 


2vr J 


2?r 



(2-84) J Zr, /K^)l(0, -^)dv ld v 2 , 

2-7T 

and it can be seen to be asymptotic to 

/ f(v U V 2 ) - -— T- .. 2 + 



2vr J-ooJ-oo ' V i( 27r ) 3,y l u 2 i(27r) 3 u|(wi - v 2 ) 

& —2m{v\—Wi) g— 27ri(«i— Ua) \ 

(2.85) + - , 3 ? y2 dv!dv2. 

l^Zir^VilVi — V2) l{2TT) a V 2 {V\ — V2) J 

All of the other terms containing a variation of I(cti, a 2 ; /3) can be obtained from (|2.83p and (|2.85p 
by a simple swapping of v\ and v 2 or by changing the sign of these variables. 

After combining like terms, the exponential terms in the integrands of (|2,83p and (|2.85p and the 
other similar integrals sum to 

g— 2nivi g27rii>i ^2ttivi ^—2isiv\ 



+ .. 5 o, r~ .. 5 o + 



4i7r 3 -i; 2 (?; 2 — 4i7r 3 ?;f (v 2 — ui) 4iir 3 v 2 V2 Mt: 3 v\v2 

e -27r«t) 2 g27rit)2 g -2-7™ 2 e 27riD 2 



+ 



4i7T 3 f|(wi — U 2 ) 4f7T 3 f|(t;i — t> 2 ) Ai7T 3 VlV 2 4iTT 3 VlV2 

e 2m(v\—V2) e -2m(v\-V2) g — 27ri(t;i — «a) g 27ri('Ui — W2) 

+ 5 ; 7? - 5 , 777 + 



4i7r 3 T;i('i;i — V2) 2 4iir 3 vi(vi — V2) 2 4i7r 3 u 2 (t>i — t> 2 ) 2 4i7r 3 i> 2 (vi — t> 2 ) 2 
sin(27r?Ji) sin(27rvi) sin(27r?; 2 ) 



27r 3 uf(v 2 — i>i) 2tt 3 v 2 V2 2tt 3 v 2 {v\ — v 2 ) 

^ 2 86 ^ sin(27rt> 2 ) sin (2tt(^i - t> 2 )) ^ sin (2tt(vi - v 2 )) 

2tt 3 viv 2 2-k 3 v\{v\ — V2) 2 27r 3 u 2 (vi — V2) 2 

Taking a common denominator, we arrive at 

- sin(27rui) + sin(27ru 2 ) + sin (2tt(vi - u 2 )) 
(2-87J — -3 -; r , 

2lT V\V2\V\ — V 2 ) 

and use the identity 

(2.88) sinx siny sin(x — y) = -(sin(2x — 2y) — sin(2x) + sin(2y)) 
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to show that all the terms in fl2.78j) containing I(a\, ot 2 ] /3) sum to 

T T f°° f°° , , sin (irvi) sin (70*2) sin U(vi - v 2 )) 
2ir 2itJ_ 00 J_ 00 irvx irv 2 ir(vi - v 2 ) 

The I\{a; /3) terms are simpler. We start with 

TL TL 

(2-90) -1- / 2 ; L [ 2 ^f(v 1 ,v 2 )h(0;^)dv 1 dv 2 . 

Picking out just the relevant terms from (12. 3p . the leading-order contribution to (12.90P is contained 



m 



TL TL 

1 fW rw f^£^M)'^ 



L 2 2ir l_TL J_TL 

2tt 2tt 



(2.91) +(^Ly~Q(i + 2np)C(l - ^p)A(^p) W d Vl dv 2 . 

Expanding the zeta functions around their pole and keeping just the leading-order term, we then 
perform the integration over T and arrive at the asymptotic result for (|2.90p for large T 

rp p /-OO f-CO / ^ g — 27rj«J2 \ 



We combine this with 

TL TL 

^ /*°° /*°° / \ ^2wiv2 

( 2 - 93 ) ~-^ lo g^/ / /( W i^2)(7^r^-7^7T2 )^i^2, 

and the sum gives 



2tt s ^2l^J_ 00 J_ 00 J ^ ^, Z/ V(2^2) 2 (2^ 2 ) 2 
T T [°° [°° f ! cos(2^2)\ 

Since Ii(a;/3) only depends on the sum of its two arguments, we see immediately that in the 

sin 2 [wivi— V2)) sin 2 (7ri)i) sin 2 (71-^2) 1 

T — > 00 limit the ii terms are responsible for W 2(vi-V2) i iPv 1 -rflv 2 m fl 2 - ' ' P - 



3. Random Matrix Theory 



We now use a very similar method to that in Section [2] to derive the triple correlation of eigen- 
values of random unitary matrices. Of course, there are more elegant methods to do this in random 
matrix theory (see Section I3.1|) but the point of Section 13.21 is that it helps to illuminate the 
preceding calculation of the triple correlation of the Riemann zeros. 
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If X is an N x N matrix with complex entries X = (xjk), we let X* be its conjugate transpose, 
i.e. X* = (yjk) where yjk = XfcJ. X is said to be unitary if XX* = I. We let U(N) denote the 
group of all N x N unitary matrices. This is a compact Lie group and has a Haar measure which 
allows us to do analysis. 

All of the eigenvalues of X G U (N) have absolute value 1; we write them as 
(3.1) e i9l ,e i92 ,...,e idN . 



For any sequence of N points on the unit circle there are matrices in U(N) with these points as 
eigenvalues. The collection of all matrices with the same set of eigenvalues constitutes a conjugacy 
class in U(N). Thus, the set of conjugacy classes can identified with the collection of sequences of 
N points on the unit circle. 

Weyl's formula asserts that for a function / : U(N) — > C which is constant on conjugacy classes, 
(3.2) I /(X)cfflaar=/ f(6 1 ,...,9 N )dX N , 

JU{N) J[0,2n] N 

where 

1 2 d9\ . . . d0N 



\e ie - 



7V!(2tt) 



N 



(3.3) dX N = \{ 

Since N will be fixed in this paper, we will usually write dX in place of dX^. 
The characteristic polynomial of a matrix X is denoted Ax(s) and is defined by 

N 

(3.4) A x (s) = det(J - sX*) = JJ (1 - se'^). 

71=1 

The roots of Ax(s) are the eigenvalues of X. The characteristic polynomial satisfies the functional 
equation 

N N 



Ax(s) = (s) N l[e- i9 "l[(l-e ie "/s) 



n=l n=l 

(3.5) = (-1) N detX* s N A x *(l/s). 
Note that 

(3.6) 8 %L {8) + l£l(l) =N . 

Ax s Ax* s 

These characteristic polynomials have value distributions similar to that of the Riemann zeta- 
function and form the basis of random matrix models which predict behavior for the Riemann 
zeta-function based on what can be proven about A. Some care has to be taken in making these 
comparisons because we are used to thinking about the zeta-function in a half-plane whereas the 
characteristic polynomials are naturally studied on a circle. The translation is that the 1/2-line 
corresponds to the unit circle; the half-plane to the right of the 1/2-line corresponds to the inside 
of the unit circle. Note that Ax(0) = 1 is the analogue of linv^oo ((a + it) = 1. 
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We let 

(3-7) z(x) = -J-j. 

1 — e x 

In our formulas for averages of characteristic polynomials the function z(x) plays the role for random 
matrix theory that £(1 + x) plays in the theory of moments of the Riemann zeta- function. 



3.1. Triple correlation by Gaudin's Lemma. Let f(x,y,z) be a smooth function which is 
periodic with period 2ir in each variable. We want a formula for 

(3.8) T 3 (/):= f f(8 n ,e j2 ,e j3 )dX, 

where the sum is for distinct indices ji, j 2 , js- It is a standard result in random matrix theory that 
by Gaudin's Lemma (see, for example, [21] or [9]), we have 

(3.9) T 3 (/) = f f(6i, 9 2 ,6 3 ) det S N (6 k - 9-) d9 3 d6 2 d6 l 



where 



(3.10) S N (0) 



sin^ 
sinf 



3.2. Triple correlation via the ratios theorem. We now produce an alternate method of 
calculation of the triple correlation for eigenvalues of unitary matrices. This method mirrors that 
produced earlier in the paper for the Riemann zeros, but many steps are cleaner and more obvious in 
the random matrix case, not to mention the fact that they are all rigorous. Therefore the calculation 
in this section serves to clarify and support the previous calculation of the triple correlation of the 
Riemann zeros. 

Let 

N 

(3.11) g(z) := A x (e z ) = ]J (l - e z e~^) . 

3=1 

Then, since g(z) has zeros at z = i9j + 27rim, m € Z, by Cauchy's theorem we have for an arbitrary 
holomorphic, 27ri periodic function /, 

(3.12) = jL j^M/(./fl d, = ± f/^mii) dz, 

where C is a positively oriented contour which encloses a subinterval of the imaginary axis of 
length 2tt. We choose a specific path C to be the positively oriented rectangle that has vertices 
5 — 7Tz, 5 + iri, —5 + ni, —5 — iri where 5 is a small positive number. Note that, by periodicity, the 
integrals on the horizontal segments cancel each other. Applying this three times, and replacing 
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each variable by its negative, we have (using the fact that C is unchanged when z — > — z), 

N N N 



(3.13) M 3 (f) := [ ^^/(Oi^i.y^ 

Ju(N) n=lj2=lj:j=1 

Ju(N) Ax Ax 



(2vri) s 




CJCJC JU(N) 



xf(izi,iz2,izs) dz% dz2 dz\ 

for a three variable holomorphic periodic function /. Notice that M%{f) is like T^(f) except that 
it is a sum over all triples of indices between 1 and N instead of over distinct indices. 

Let C- denote the path along the left side of C from —5 + ni down to —5 — tri and let C+ denote 
the path along the right side of C from 5 — tri up to 5 + ni. Ignoring the integrals over the horizontal 
paths (because their contribution is 0) we take each variable Zj to be on one or the other of the two 
vertical paths C_ or C + . In this way our expression can be written as a sum of eight terms, each 
term being a triple integral with each integral on a vertical line segment either C_ or C+. These 
eight integrals are analogous to J±, . . . , Js in Section T2.31 As we did for the Riemann zeta function, 

for each variable Zj which is on C_ we use the functional equation (|3.6p to replace e~ z i (e~ z i ) by 
A' 

N — e z i -f^(e Zj ). In this way we find (using X -1 = X*) that 

( 3 - 14 ) = 7^3 E / / / L,M^ N+ ^ Z] iB^^) dx 



3=1,2,3 

xf(iz\,iz2,iz3) dz3 dz2 dz 



We next compute the averages over X £ U(N) which appear in the above equation. 

We need two instances of the ratios theorem (see |12| or [13] for statements and proofs of this 
theorem). Recall the definition of z(x) from (|3.7p . 

(A) Let m > 0. 

(3 15) f Ax(e- a )^x*(e-P) dx = z (a + P)z( 1 + 5) ^ N{a+/3) z(-(3 - a)z( 7 + 5) 
J U(N) A x (e-T)A x .(e- 5 ) z(a + 6)z(J3 + 7) + S)z(-a + 7) ' 

(B) Let St^S^SM > 0. Then 

kx(e-^)kx(e- a *)hx*(e-P) dx _ z(a x + (3)z(a 2 + (3)z( 7l + 5)z( l2 + 5) 
U(N) A x (e-^)A x {e-^)A x *(e- s ) z{a x + 5)z{a 2 + S)z(/3 + lx )z{fi + 72) 

+e _7v (ai+/ 3) - ai)^(a 2 ~ «i)^(7i + S)z(-y 2 + g) 

z(— fi + <5)z(a 2 + <5)z(-ai + ji)z(-a\ + 72) 

(3 16) I c - N(a2+0) Z<K ~^ - a ^ z ^ ai ~ " 2 ) z fa + ^^(^ + 8) 

z(—(3 + 6)z(ai + S)z(— ct2 + 7/1) z(— a 2 + 72)' 

For application to triple correlation, we need averages of logarithmic derivatives of the charac- 
teristic polynomials. Differentiating the above formulas leads to the following two identities, which 
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are the random matrix version of (|2.3f) and (12. 2p respectively, remembering that z(x) plays the role 
of C(l + x) and N the role of log The products and sums of primes in (12. 3ft and (|2.2I) do not 
have any counterpart in random matrix theory. 

Proposition 3.1. If$ta,$t./3 > 0, then 



J(a;/3): = e — " / ^(O^e^X 
JU(N) A X Ax* 

(3.17) = (£\ (a + P) + e- N( - a+ ^z(a + P)z(-a- 0). 

Proposition 3.2. Let 9?a 2 , 5R/3 > 0. Then 

J( ai ,a 2 ;P): = - e "— * f ^k( e — )^( e — )^(e^)dX 

JU(N) A X A X Ax* 

= e- Ar(Ql+ «z(ai + /3)z(-ai - /?) ( -(a 2 - ai) - -(a 2 + /?) 



(3.18) +e -m*2+P) z ( a2 + /J) z (_a 2 - p) ( -(at - a 2 ) - ~(ai + (5) 

\z z 

Note that the right-hand-side has a simple pole when ct\ = —(3 and when a 2 = —/3 but is analytic 
when ai = a 2 . We also remark that non-constant integrals with no X, or no X*, are 0; for example, 
if Ka, > 0, then 

(3.19) / £*(0 = / ^l(e-) d* = / ^(e' Q )^(e-^ rfX = 0. 
■/[/(AT) A x JU(N) A X* ./[/(AO A x Ax 

Thus, letting dG be a shorthand for f(i Zl ,iz 2 , iz 3 ) dz 3 dz 2 dzi, we have, upon expanding (|3. 14[) . 

(3.20) -(27Ti) 3 M 3 (f) = N 3 [ [ [ dG+ [ [ [ J( Zl ,z 2 ;-z 3 ) dG 

JC- JC- JC- JC+ JC+ JC- 

+ J(z 2 ,z 3 ;- Zl ) dG+ / / J(z 1 ,z 3 ;-z 2 ) dG 

JC- Jc+ JC+ JC+ JC- Jc+ 

{J{- Zl , -z 2 ; z 3 ) + N(J(- Zl ; z 3 ) + J(-z 2 ; z 3 )) dG 
{J{-z 1 ,-z 3 -z 2 )+N{J{-z 1 -z 2 ) + J{-z 3 -z 2 )) dG 
{J{-z 2 , -z 3 ; Zl ) + N(J(-z 2 ; Zl ) + J(-z 3 ; Zl )) dG. 



C- JC- JC 



C- JC+ JC 



C+ JC- JC- 

This is very similar to (|2.53[) , except that we have retained three variables instead of working with 
the differences z 2 — Z\ and z 3 — Z\ . 
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Now, by the inclusion-exclusion principle, 

(3.21) T 3 (f) = M 3 (/)- / > {f{p h ) + f(0n 

, 6j 2 , 0^ ) +/(%!, 0j 2 , #j 2 ) ) cLX" 

+2 / £/(0^,%)dX 

The pair-correlation sums are evaluated much as above. For example, 

(3.22) f J2f(6 jl ,6 h ,6 j3 )dX= 1 (- f [ J(-zi\ z 3 )f(iz 1 ,iz 3 ,iz 3 ) dz 3 dzi 



- / J(zi;-z 3 )f(izi,iz 3 ,iz 3 )dz 3 dzi 
Jc+ Jc- 

+N 2 J j f(izi,iz 3 ,iz 3 )dz 2 dzi^. 

The one-correlation sum is 

(3.23) / V f{0j,0j,0j) dX = ^r [ f(iz, iz, iz) dz. 

JU(N) 77ir 27rz JC- 



Now we move all of the paths of integration over to the imaginary axis. When we do this we 
encounter some poles and have to use principal value integrals. For example, 



(3.24) III J^i,z 2 ;-z 3 ) dG 

/ni p r 
/ / J(zi,z 2 ; -z 3 ) dG 
-iii JC + JC- 

/ni p p 

lim I I J(zi,z 2 ; -z 3 ) dG 

z 2- z l l> e 

/wi P P 

lim / lim / i-™,™] J(z±, z 2 ; —z 3 ) dG 

. %i 6^0+ J 6^0+ J \z 3 -z 2 \>5 

z 2- z ll> £ \z 3 -z!\>6 

/7ri p 
lim / Res ( J (zi,z 2 ;—z 3 )f(izi, iz 2 ,iz 3 )) dz 2 dzi 

-TTi <^0+ J [-Ti.TTi] 2 3 =Z1 
z 2- z ll> E 

/7T2 P 
lim / Res (J(zi,z 2 ;—z 3 )f (izi,iz 2 ,iz 3 )) dz 2 dzi. 

-TTi £-0+ J , I-™,**] 23=^2 
z 2~ z ll> e 

The first equality follows because there are no singularities when we move z\ onto the imaginary 
axis. In the second, there are again no singularities when we move z 2 onto the imaginary axis, 
but anticipating what comes next we choose to write this integral with a limit as e — > + . Finally, 
the third equality takes into account the poles at z\ and z 2 when we move z 3 onto the imaginary 
axis. To determine the signs of the three terms in the final equality, remember the the contour C_ 
is oriented downwards and just to the left of the imaginary axis. Thus when it is moved onto the 
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vertical axis it wraps the singularities in the positive (anti-clockwise) direction. Also the residues 
are multiplied by vri instead of 2m because we moved the path onto a path that goes right through 
the singularities, i.e. we have not crossed the poles but rather moved on top of them, and so we 
have only one-half the usual contribution of a residual term and the integral in z 3 that remains in 
the first term of the final equality is a principal value integral. 

Now it is easily calculated that 



Res J(zi,z 2 ;-z 3 ) = -J{z 2 ;-z\) and Res J(z±, z 2 ; — z 3 ) = —J(zi;— z 2 ) 



Thus, we have 



(3.25) III J(zuZ2-,-z 3 ) dG 

J C+ J C+ J C— 

/iri p p 

lim / lim / i-™,™] J(zi, z 2 ; — z 3 ) dG 

_ ni e-*0+ J 8^0+ J \z 3 -z 2 \>6 

z 2- z l!> £ \z z - zi \>6 

/Tri r 
lim / J(z 2 ; — z\)f(iz\, iz 2 , iz\) dz 2 dz\ 

-Tfi e - 0+ J , [-^'^ 
z 2- z ll> e 

/TTi P 
lim / J(zi;—z 2 )f(izi,iz 2 ,iz 2 )dz 2 dzi. 



-m ' - - 1 

z 2- z ll> e 



A slightly different calculation gives 



(3.26) SIS J(z^z 2 -- Zl ) dG 

J C— J C+ J C-\- 

/7VI P P 
/ / J(z 3 ,z 2 ; -zi) dG 
-iri Jc+ Jc + 

lim / / J(z 3 ,z 2 ; -zi) dG 

-TTi «-0+ J J C 

\ z 2- z l\> e 

+iri / Res ( / J(z 3 , z 2 ; -zi) f(izi, iz 2 , iz 3 ) dz 3 ) dz\ 

J -TTi 2:2=21 V JC+ J 



/wi P P 

lim / / J(z 3 ,z 2 ; -zi) dG 

.jri e-0+ V J c 

\ z 2- z l\>< ! 

/ni p 
J J(z 3 ;-z 1 ) f(izi,iz 1 ,iz 3 ) dz 3 dz\ 
-TTi J Cj- 
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after moving z\ and z 2 onto the imaginary axis. Now we move z 3 onto the imaginary axis and 
obtain a residual term from each of the two integrals above. Thus, 

(3.27) / / f J(z 3 ,z 2 ;-z 1 ) dG 

JC- Jc+ Jc+ 

/7ri p p 

lim / lim [-^i] J (23, z 2 ; — z\) dG 

-iri e— 0+ J <5-»0+ J |* 3 -* a l>* 



-7ri e — >yJ J 1 ""'""J o — >\j 1 j 123- 

l z 2~ z ll> E \z 3 - zi \>5 



TV I 



-ni / lim / Res (J(z 3 , z 2 ; — z±)f(izi, iz 2 ,iz 3 )) dz 2 dz\ 

e^0+ J [-"«,«] 23=21 
l z 2~ z ll> E 

(•TTl r 

-iri I lim / J(z 3 ; — z\) f(izi,iz\,iz 3 ) dz 3 dz\ 

->0+ J [-"'>"•] 

l z 3- z ll> e 



— TTl 



(iri) 2 I Res [J(z 3 ;-zi)f(iz\,izi,iz 3 )) dz 



Now 



Res J(z 3 , z 2 ; — z\) = —J(z 2 ;—zi) and Res J(z 3 ;—z\) = —N. 

23=21 2 3 = Zl 



Thus, we end up with 



(3.28) / / / J(z 3 ,z 2 ;- Zl ) dG 

J C— J C_)_ J c+ 



/7T1 r r 

lim / lim / i-™,™] J(z 3 , z 2 ; —z\) dG 

-Trie— 0+ J <5^o+ 7 l»3-*al>* 



l z 2~ z ll> E |z 3 - Zl |>i 



77/ 



-tt« / hm / J(z 2 ; — z\)f(izi, iz 2 , iz\) dz 2 dz\ 

j e~^0+ J [-*■»,*■>] 
l z 2- z ll> E 



77/ 



-7ri / lim / J(z 3 ;—zi) f(izi,izi,iz 3 ) dz 3 dz\ 

l z 3~ z ll> E 



7T/ 



— (iti) 2 N I f(iz,iz,iz) dz. 
By a change of variables we have 

(3.29) ( [ [ J(-z 3 ,-z 2 ; Zl ) dG = - [ [ [ J(z 3 , z 2 ; - Zl ) dG. 

JC+ JC- JC- JC- JC+ JC+ 

We can also make use of the symmetry J[a\, a 2 ; (3) = J(a 2 , a±; P). In this way we can use f|3.25j) 
and (|3.28j) to replace all of the integrals involving a three-variable J in terms of principal value 
integrals. 
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In a similar way we rewrite, for example, the two variable J- integral as 

(3.30) / / J(zy,-z 3 )f(izi,iz 3 ,iz 3 )dz 3 dzi 
JC+ JC- 

lim / J(zi;-z 3 )f(izi,iz 3 ,iz 3 )dz 3 dzi 

r j e^0+ J 

l z 3~ z ll> e 

/Tri 
f(iz, iz, iz) dz. 

Using (|3.14p - (|3.30p . after the substitutions izj = 0j, j = 1, 2, 3, and some simplification, we have 
WHtAs / lim ,/, , lim J (j(i0 1 ,i6 2 ;- l e 3 ) + J(ie 1 ,ie 3 ;-i9 2 ) 

+j(ie 2 ,ie 3 , -io l ) + j(-ie u -ie 2 -ie 3 ) + J(-ie u -io 3 - w 2 ) + J{-io 2 , -io 3 \ ^i) 

+N{j{-iO x ;i0 3 ) + J{-i6 2 -i9 3 ) + J{-i9r,i9 2 ) + J{-i6 3 -i0 2 ) 

(3.31) + j(_ i 2;i 1 ) + J(-ie 3 ;i0x)) +N 3 ^f(6 1 ,6 2 ,8 3 ) dB x d0 2 d9 3 , 

all of the pair-correlation and one-correlation terms having cancelled. The principal value integrals 
with the limits in e and 5 are no longer needed because it can be checked, using Mathematica for 
example, that the integrand here is 

(3.32) = det S N (9 k - ^)/(0i, 2 , #3) 



which is entire. Now our formula agrees with (13.91) . 

This concludes our alternate proof of triple correlation for unitary matrices. 
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